THE MAXIMAL VELOCITY OF A PHOTON 



AVY SOFFER 

Abstract. 

We estimate the probability of a photon to move faster than light, under the 
Einstein dynamics which, unlike the wave equation or Maxwell wave dynamics, has 
singular dispersion relation at zero momentum. We show that this probability goes 
to zero with time, using propagation estimates suitably multiscaled to control the 
contribution of low frequency photons. 



Section 0. Introduction 

Let F(X > c) stand for a smoothed out characteristic function of the interval 
A > c. We say that a quantum system, with a configuration coordinates x, obeys 
a maximal velocity bound, if for its any evolution i(j(x,t), localized in a compact 
energy interval , there is a positive constant c s.t. 



>c)1>[x,t)\\ La ->• 0, 



\x\ 

T 

as \t\ — > oo. It is fairly easy to prove such bounds for non-relativistic iV— particle 
quantum systems (see [SS, HSS]). However, such bounds are still open for photons 
interacting with such systems. The goal of the present paper is to prove such 
bounds for a model of a single photon with an effective interaction described by a 
potential V(x). This model is given by the time-dependent Schrodinger equation 



dtp 

(0.1a) i— = Hi/j VU=o = ^o, 

with the Schrodinger operator H on L 2 (IR 3 ), given by 
(0.1b) H := (-A) 1/2 + V(x). 
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This Hamiltonian is a simplified, scalar field version of the restriction of the photon 
field in QED, to the one particle sector. We think of the equation (0.1) as a 
model describing dynamics of a single photon. Moreover, it is used as a laboratory 
for developing the methods to to study propagation with singular dispersion at 
zero energy, needed in tackling the full non-relativistic QED problem. Notice that 
the equation, without the interaction term, is relativistic invariant. Moreover, by 
multiplying the equation by — idt — \p\, we get the wave equation. However, since 
the initial data for the derived wave equation, is not localized: 



we can not use the finite propagation property of the wave equation for this model. 
We make the following assumptions: 

(i) V(x), x • VV(x) are sufficiently regular and decay faster than 0(|x| -2 ) at 
infinity; 

(ii) H has no zero energy resonances or zero energy eigenvalues. 
The main result of this paper is the following theorem: 

Theorem. Let ip(x,t) be the solution of the Schrddinger equation (0.1), and as- 
sume, furthermore, that H satisfies the conditions (i) and (ii). If the initial data 
is localized in energy in some compact interval, and is such that < x > 1+e Vo has 
a bounded L 2 norm, we have the following maximal velocity bound: 



for R^> 1 and \t\ — > oo. 

As was mentioned above, the proof of such a propagation estimates for general 
N— body Hamiltonians is quite elementary, see e.g. [Sig-Sof, HSS]. See also [BFS 
and cited references]. However, the fact that the one particle Hamiltonian of a 
massless photon is singular at zero energy, prevents a direct extension of the general 
theory of propagation estimates to the field, see however [BFSS, FGS]. 

The proof is based on the construction of propagation observables, on each 
(dyadic) scale of the energy. Suppose that the length of the momentum, p = — iV, 
is localized in [2 _n_1 , 2 _n ], n > 0. Consider then the operator 



^(o) = MV(o), 



(0.2) 




F(A/Rt2~ n > 1) 



defined through the spectral theorem. 



(0.4) 



A = - 



i/2(x- V + V -x), 
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and R> 1. 

Then, we have the basic identity, derived from the Schrodinger equation: 
(0.5) d t < 1>(t), Fi(>(t) >=< m, {*[H, F] + d t F}if;(t) > . 

Here, F stands for any self-adjoint operator, for which the differentiation can be 
justified, for the chosen ip(t). To get a useful propagation estimate, we use the 
idea of negative propagation observables, [Sig-Sof ] : Suppose F stands for a family 
of time dependent operators, bounded from above (e.g., negative), and such that its 
Heisenberg derivative, defined above, is positive, up to integrable (in time) correc- 
tions. We will then obtain, upon integration over time, the following propagation 
estimate: 

(0-6) jf ||^(t)|| 2 ^<C||^(0)|| 2 . 

Here, we use that the Heisenberg derivative, D H F, is given by a positive operator, 
plus integrable corrections: 

(0.7) D H F = i[H, F] + d t F = B*B + R(t). 

In general, we need to use phase space operators which, while not pseudo-differential,( 
will have suitable phase space support. If we choose for F, the function of A, de- 
fined above, as our propagation observable, we need to estimate the commutator 
of H with this F, suitably localized in some energy shell. One can then use the 
following basic commutator expansion Lemma: 

(0.8) i[H, F(A)} = F'(A)i[H, A] + R([[H, A], A]), 

where the remainder is given by an explicit multiple integral, involving the group 
generated by A, and, the double commutator above, of H with A. By the above 
expansion, we obtain, through symmetrization, that the above propagation observ- 
able, has Heisenberg derivative, that is smaller than: 

(0.9) A G l {A /t) + R(t,n). 

Here, G n stands for a bump function of fif ^_„ around 1. We will get a useful 
estimate, if we can control the remainder term, R(t,n), by an integrable function 
of t, which is also well behaved as n tends to infinity. The leading term in the 
expression (0.9), comes from the time derivative of F term, of the Heisenberg 
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derivative. The first term in the Heisenberg derivative, the commutator, has an 
expansion, beginning with: 

on 

(0.10) F\A/t)i[H,A]-, 

which can be shown to be much smaller than the leading term, by using the en- 
ergy localization around 2 _n . To see that, we use that the commutator i[H, A] is 
bounded by c\p\, for some finite constant c, by the use of the uncertainty principle 
and our decay and regularity assumptions on V. Then, we need to show that energy 
localization implies a similar momentum localization. This can not be done using 
standard localization arguments, since the derivative of the localization function 
grows like 2 n ! To this end, a completely different argument is used: One proves 
that under generic spectral assumptions on H, HP C (H), dominates a constant times 
P c (H)\p\P c (H). P C {H) stands for the spectral projection of H on its continuous 
spectral part. 

It then follows that: 

(0.11) E n (H)\p\E n (H) < c2- n E n (H) 2 . 

Finally, we need to show that the remainder R(t,n), higher commutator terms in 
the expansion of the Heisenberg derivative, are integrable in t. This is the most 
involved step, without energy localization the formal expression for the remainder 
is given by a divergent integral. The way to estimate this last remainder term, is to 
use the fact, that the group generated by dilations moves the support of the energy 
localization functions. Consequently, the integrals over the group actions is limited 
to finite domains. The remainder term comes from the following expansion: 

E n i[\p\, <S> n ]E n = E n i[\p\ V 2 , ^ n ]\p\ 1 / 2 E n + Enlp^HM 1 / 2 , <S> n ]E n 

(0.12) = E n \p\ 1/2 K(A/t)\p\ 1/2 E n ^- t + 2m n \p\ 1/2 R2(A/t)E n , 

with 

A 

$ n = Ei F( > l)Ei 

and with R 2 (A/t) given by 

R 2 (A/t) = \ [ d\F n (\y XA / Rt [ X ds f du e~ luA / Rt hp] 1 / 2 e iuA / Rt (Rt)- 2 
2 J Jo Jo 2 

(0.13) =- A j ' d\F n {X)e lXA / Rt J\s J S du e- u l Rt \p\ x l 2 (Rt)- 2 . 
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If we try to bound the expression for R 2 , by taking the norm of the first term 
on the rhs, as it was done in past works, we lose a factor of 2~ n / 2 , coming from 
localizing the |p| 1//2 factor. On the other hand, we can not directly estimate the last 
expression on the rhs of equation (0.13), since the integrand grows exponentially, 
while F n decays slower than exponential, being the Fourier transform of a compactly 
supported function. To this end, we use the fact that the dilation group, generated 
by A, changes the support of functions of \p\, or H: 

(0-14) E In (\p\)e iXA E In (\p\) = 0, 

for |A| > In 2. Then, we use the mutual domination of \p\, H: 

(0.15) P C (H)H < cP c (H)\p\P c (H) < dP c (H)H, 

for some positive constants c, d. 

Combining equations (0.14), (0.15), we can then show that the integration on A 
is limited to a compact domain, in equation (0.13). Collecting all of the above, we 
get estimates of the form: 

(0.16) E n (H)^G 2 n (A/t)E n (H) e L\dt). 

This estimate is then jacked up by the use of the propagation observable ^F n {A/t), 
to obtain, 

(0.17) E n (H)F n (A/t)E n (H)j e L\dt). 

In the next step of the proof, we estimate, using the above, the following operator: 

E n (H)F(^ > R)E n {H). 

We write, F(M > R) = F(M > R)F n (A/t) + F(M > R)F n (A/t). 

F n + F n = 1. 

The first term of the above decomposition, goes to zero, as time goes to infinity, by 
the above propagation estimates, on F n (A/t). 

So, we need to show that the F n term also goes to zero. 

This is formally true , since it consists of a product of two operators, with disjoint 
classical phase-space support, on the energy shell 2 _n . 

Again, the proof of this property necessitates the use of new phase space lo- 
calization arguments. In the final step of the proof, we sum over all n, and in the 
process we lose some powers of 2~ n . These are compensated by requiring the initial 
data to be localized in x, and by using that(up to 2) negative powers of \p\, are 
bounded, up to a constant, by positive powers of \x\. 
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Section 1. Propagation Estimates 

Our goal in this section is to prove the following key propagation estimate, as 
sketched in the introduction: 

\\F n (A/t)E n mf T <c\\E n m\\ 2 - 

and with 



(1-lb) Fn (A/t) = F( 1 ^ =K >l). 

E n stands for the operator E In {H). We use propagation estimates, with the prop- 
agation observables 

(1.2a) $ n = E In (H)F > A Eln (H), 



(l-2b) $ n = E In (H)jF > lj E In (H), 

where, n = 0, 1, • • -oo, I n stands for the interval [2~ n_1 , 2 _n ], E In is the charac- 
teristic function of I n and A is the dilation generator, 

(1.3) A = ^(x-p + p-x). 

Note that ^ > R on the support of $ n , and R < ^ < 2R on the support of 
E n F'E n . Furthermore, 
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The main propagation estimate is based on showing that: 



Ej i 



\p\ + V,F 



A 



Rt2~ 



> 1 



(1.5a) 



= E n {i\p\ 1 ' 2 [M 1/2 ,F] +i [M 1/2 ,F] M 1/2 + i[V, F]} E n 

= E n {bl^^H^ + |p| 1/2 |^Q 2 (^A)bl 1/2 } E n + or 2 ), 



where F 2 ~ j^F' and, E n Q 2 {A/t)E n is of order 1 for t > 2 n /R, and of order 
2~ n , for t < 2 n /R. Then, we show that 

(1.5b) E n \p\ 1/2 F 2 \p\ 1/2 E n < (c/Rt)E n G 2 n E n , 

with 



Gl = F{ 



Rt2~ 



1). 



Letting F n (A) = F(2 n A > 1) and using that R2~ n F^R > 1, we 

obtain 

J d t <iP u E In (H)F n (^- t )E In (H)ip t >dt= J <MH^n]ipt>dt 
(1.6) +2 n R~ 1 £ < ^ u E In {H){^)F' n E In {H)^ t > dt. 

Assuming that (1.5) holds, using (1.6) we can then prove: 

Theorem 1.1. Under the previous assumptions on H , and assume that (1.5) holds, 
we have the following propagation estimates: 



(1.7a) 



\K(A/t)E n m\\ 2 j < c(R)\\E n m\\' 



T dt 

\\F n (A/t)E n m\\ 2 -r 
\ 1 



(1.7b) 



+ < ip(T), E n —F n (A/T)E n ip(T) > < c(R)\\ < A >^ 2 £^(0)|| 2 
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proof. 

To prove (1.7a), we observe, that by the estimate (1.5), combined with Propo- 
sition^^), each power of \p\ contributes a factor of 2 _n , the F term is bounded 
by 

(c/Rt)F f , which is dominated by the integrand of the second term in equation 
(1.6), for R large enough. The Q 2 term contribution to the commutator in equation 
(1.6), is bounded by our assumptions , following estimate (1.5): 

||£«bl 1/2 || < c2" n/2 , 

by Proposition (2. 4), which eliminates one power of 2 n ; then, integrating 1/t 2 over 
the region of t > 2 n /R, eliminates a power of 2 n . The integral over t < 2 n /i?, is 
also bounded, since Q2 is assumed to be smaller than c2~ n . To prove (1.7b), we use 
the propagation observable (1.3b). The proof is similar to the first case: The time 
derivative term, which contributes the second term on the rhs of (1.6), is , as before 
the dominant term. The Q2 term is controlled as before (the only difference is that 
now F is replaced by d\F(X), in the estimates on Q 2 . As for the commutator term, 
we have now 

1 A 2 n 

(1.8) i[H, $ n ] = E n -i[H, A]F n (A/t)E n + E n —F^(A/t)i[H, A]—E n + E n Q 2 E n . 

The Q2 term is estimated as before. The second term on the rhs of 91.8) is also 
estimated as before, since the factor A/t ~ 2~ n j R, on support F' . However, the 
first term is unbounded, so the bound for t < 2 n /R, Proposition (C), does not 
apply. Instead, we break the F to two parts, one decaying at infinityand sharply 
localized in Fourier space, and one constant at infinity. The first part is estimated 
as before, and the second part is directly controlled by elementary size estimates, 
see Proposition (D).B 

There are new difficulties in completing this argument, compared with the usual 
case, without dyadic energy localization. 

First, we need to minimize the number of powers of 2 n , coming from expand- 
ing the function F n . Then, we need to trade positive powers of the momentum( 
derivative operator) p = — zV, for powers of 2~ n . 

Finally, to control the remainder term in the Commutator Expansion Lemma, 
the Qi term, (or R2 term), we need to commute the derivative through the dilation 
group, which produces exponentially large factors. 

The way out of these problems involves the following arguments. To limit the 
integrations in the remainder term R2, we notice that, the dilation group moves the 
dyadic energy interval, away from its original support. Hence, for large enough value 
of the group parameter, A, the fact that our propagation observable is localized 
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on the dyadic interval, from both sides, gives an extra decay, that cancels the 
exponential growth factor. This is shown in detail in the subsection "The term 
R 2 n . 

To get the 2~ n factor from the momentum p, we prove some propositions about 
the properties of the operator H, which might be of independent interest, (see Propo- 
sition (2.4)) Specifically, we show, that in three (and higher) dimensions, if there 
are no zero energy resonances and eigenvalues, then H and \p\ dominate each other, 
up to a multiplicative constant, on the continuous spectral subspace of H. 

These estimates are the key to getting the right minimal powers of 2 n , from the 
various propagation estimates and phase space localizations. 

To prove the necessary estimates, to fill in the details of the above theorem, and 
the assumptions on Q 2 , we estimate the Heisenberg derivative of the observables 
defined in (1.3), on the time intervals t > 2 n /R, and t < 2 n /R. This is the content 
of the following four Propositions A-D. 

We begin by deriving and estimating the R 2 term: 



The R 2 (A/t) term. 



Proposition A. 

For $ n as in (1.3a), we have: 



(1.9) / < E n \p\ 1 / 2 Q 2 \p\ 1 / 2 E n >dt< --||£ n V>(0)|p 

J2 n /R K 



Proof. 

Direct application of the commutator expansion lemma gives: 



p\ 1 ' 2 ,F n {A/Rt)}= [ d \F n (\y XA / Rt ±- [ e-^^HUlpMe^^'ds 
i J Rt J L J 

(1-lOa) = iF^A/t^p^ + RM/tM 1 ^ = iF^A/t^p^ + RM/t), 
where we used that 

(l.iob) i[A\p\ 1/2 } =p-v p \p\ 1/2 = \\p\ 1/2 - 

R 2 (A/t) = ~ I d\F n (X)e lXA / Rt t ds f du e -™A/mh p] i/ 2 ^uA/m^y 
2 J Jo Jo 2 
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(1.11) = ^J d\F n (X)e lXA / Rt £ ds jf du e- u / Rt \p\ 1 / 2 (Rt)- 2 . 

In general, this integral blows up at infinity, due to the fact that e u l Rt grows expo- 
nentially fast, while F n (X) decays faster than any polynomial, but not exponentially, 
since F n (X) is the Fourier transform of a compactly supported function. 

E n i[\p\,<f> n ]E n = £ n z[b] 1/2 , $ n ]M 1/2 £n + Kb| l7 ^[b| 1/2 , $n]E n 

(1-12) =EM 1/2 K{m\p\ 1/2 E n ^ t + E n \p\ 1 ' 2 R 2 {A/t)E w 

with R 2 (A/t) given by (1.11). 

E n \p\ l ' 2 R 2 (A/t)E n = E n \p\ 1 / 2 E n (\p\)R 2 (A/t)E n (\p\)E n 

+ E n O(2- n n)E n (\p\)\p\ 1 / 2 R 2 (A/t)E n (\p\)E n 
+ E n \p\ 1 / 2 E n (\p\)R 2 (A/t)E n (\p\)0(2- n n)E n 
(1.13) = J1 + J2 + J3. 

E n = 1 — 

In our case, 1 — E n = E(\p\ < 1) — E n . E n stands for smoothed E n function, and 
where we used proposition 2.2(d). 



J Jidt = J^ E n \\p\ l ' 2 E n (\p\)R 2 (A/t)E n (\p\)E n dt 

r E n b| 1 / 2 E n (|p|)g 2 (AA)E n (b|)b| 1 /2 i?n 

-Kb| 1 / 2 E n (|p|)Q 2 (AA)K(b|)|p| 1 /2^ 



T at 
1 m 2 

° 2 dt 



1 R 2 t 2 

2n ^ E ^\P\ 1/2 ^mQ2{A/t)E n {\p\)\p\ l / 2 E n . 

If T < c2 n , then the second term on the r.h.s of (1.14) is zero. 
Estimating Jl 

Using proposition (2.2c), it follows that the A integration (and therefore the 
other integrations) is limited to 

(1.15) |A|<ittln2. 
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Hence, 

= j^w- 2 En{H)\p\ 1,2 E n {\p\) [ dXF n {X)e iXA ' Rt 

\ m ) J\\\<Rt\n2 



(1.16) x f ds [ S due luA/m \p\ 1/2 e- tuA/Rt E n (\p\)E n (H). 

Jo Jo 

J1 = 7^ E ^P\ 1/2E n(\p\) [ F n (X)e iXA/Rt f ds f e^^du^EnQpDEniH) 

\ m ) J\\\<Rt\n2 Jo Jo 

= ^0(E n 2- n / 2 [ \X 2 F n (X)\ dX 2~ n / 2 E n ) 



\\\<Rt In 2 

' ; 0(E n 2- n (Rt) 2 E n ). 



~~ (Rt) 2 v ' n >- | 

Hence, 

f 2 n /R 1 

J Jl dt<0{-). 
If Rt > 2 n , we use instead, that 

J \X 2 F n (X)\ dX<2 2n , 

so that, 

31 = 7^2°( E ^ nE n)^ and then, 
[Rt) 



(1.17) / Jldt = 0{\, 

J2 n /R n 



/J 



The estimates of J2 and J3. 



Consider the region Rt > 2 n . The integrand to estimate, which is 
(|p| s = e ~ lAs ' Rt \p\e lAa ' Rt ) , can be written as, 

(Rt)- 2 F n (X)J(X,s,u), 
J(X,s,u) = E n {H)\p\y 2 E n \{ P \)e^ A ' R M\I 2E n m E n{H), 
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and adjoint of such term. 

First, we decompose the region of integration A to: 



IAI 



IAI 



— > m and — < m, m > In 2. 
Rt Rt 



For ^ > m we commute E n (\p\) through e lXA / Rt , to get, 

for some n(A) 7^ n , (since ^ > In 2, by assumption), and so E n E n ^x) = -E'n(A)- 
Then, 

/ ^272 / F nW / / J(X,s,u)dXdsdu 

J2 n /R R * J -00 Js=0Ju=0 



[ \F n (X)X 2 \d\-2- n / 2 -2" 



< — 2" n i?sup 

< -2- n / 2 2- 2n O(\m2 n \- k ). 
R 

In the last estimate, the factor 2~ n R comes from the t— integration. The factor 
2~n/2 comes f rom |p| 1 / 2 £' n (|p|). This factor is missing in the adjoint terms. Another 
factor, 2~ n , comes from applying proposition 2.2 to the product E^x)(\p\)E n (H). 

Finally, we note that, by construction, F n (X) is vanishing faster than any poly- 
nomial at infinity, and is approximately a constant, of order R2~ n , on an interval 
of size 2 n /R around the origin. Since |A| > mRt > m2 n , the 0{\mt\- k ) bound, 
Vfc, follows. 

Next, we consider the region ^ < m. Now, ^ < m =>• ^,^<w. 
Therefore e \s\/Rt^ e \u\/Rt < e m ^ go ^ |p| s < e m |p| . So, we pick up a factor of 
2~n f rom Ipj 1 / 2 factors, and the integration of it, s,t gives: 



Hence, 



•4<- 



A 2 F n (A)|dA < c(2 2n /i? 2 )e m 



J2 n /R 



dt 
RH 1 




J(A, s, u)F n (X)dXdsdu 



< c—2- n R2- n (2 2n jR 2 )2- 



where we gain an extra 2 n from E n E n and |p| each. 
Large |p| or if. 
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Since the momentum operator is unbounded, we need to bound it, under the 
integrals, in the definition of R2, for large values of \p\, or H. To this end, we 
replace \p\ — > Hg(H < 10) ,and then, we need to work with: 

-\p\ — >• i[A, Hg] and i[A, [A, Hg}]. 



-i[A, Hg] = (\p\ + x ■ VV)g - Hi[A, g(H)] 
= g(\p\+x-VV)-i[A,g(H)]H. 

This change is not difficult to handle as before. 

In cases we only used the boundedness of these commutators, it is straight- 
forward, since both g{H)\p\ and i[A,g{H)]H are bounded. 

In case we need to use a factor of \p\ (or If] 1 / 2 ) to get extra smallness factor 2 _n , 
it is done as before, since the g(H) factor does not change E n (H). The x ■ W term 
is better, since it is bounded by c\p\ 2 . ■ 



The region < t < 2 n /R. 



Proposition B. 

For $ n as in (1.3a), we have: 

r 2 "/ R tit 

(1.18) J \\G n (A/t)E n m\\j < c\\E n ^)f 

Here, G n , is a bump function of A/t around 2~ n /R. 

Proof. 

Now we have, 

d , / r/<I> 



where, 



m ^^r m ) = (Enm.t-'F'^A/t) E n m) 

< ~ (E n ^,2- n RF' n {A/t)E n ^) = -1 (E n iP,Fl{A/t)E n ip) , 



\F 2 (A/t)\<l and F 2 > 
F 2 (A/t) ~ 1 for - ~ R2~ n 

t 

F 2 (A/t) = for -00 R2~ n . 

t 



14 



AVY SOFFER 



i[H, $„]V(*)) = (VM, i(HE n F n E n - E n F n E n H)^(t)) 

= i (E n iP(t), (HE n F n E n - E n F n E n H)E n i;(t)) 
<2-2- n \\E n iP(t)\\ 2 . 



Hence, 



dt 



^0 



F n {A/t)E n *l){t) 
F n {A/t)E n ^{t) 



+ 2 



2 2 



£nV>(*) 



£ n V>(t) 



2" n cft 



Therefore, 

t 



Jo 



F n (A/t)E n ^(t)|| 2 + (i>{2 n /R),F n (RA/2 n M2 n /Rj) 



Improved Decay 

Proposition C. 

Under the assumptions as the above Propositions (A,B), we have the following 
propagation estimate, 

r-T 

IF 



> E n tfj(T) > < c(R)\\E n < A >V2 F n ^(0)|| 2 , 



(1.19) 

+ <^(T),£ n ^F n (^^ 
/or aZZ sufficiently large. 

Proof. We can use E n jF n (A/t)E n = $ n . The estimate of R,2(A/t) is the same 
as before, so for t > 2 n , it is done as before. The main change is that now F(X) is 
replaced by d\F(\),, and a similar proof applies. The first term in the commutator 
expansion of i[H, <fr n ], has an extra, positive term, which is however integrable over 
time by the Previous propositions (A,B), since it is supported (in phase space) on 
the support of F' n . ■ 

For < t < 2 n , since jF„ is not bounded, the proof is different. 
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Proposition D. 

r 2n / R ( a \ rit 
I "^(^> 1 )^ )II T 

(1.19) 

+ < ^(2 n /R),E n ^^F n (A > l)E n ^{t) > < c(R)\\E n < A > x l 2 K^(0)|| 2 , 

for all R sufficiently large. 

Proof. 

We write, 



where, 



and 



^4 



Fn,! = i J e- x2 d x F n (\)e-^ x d\, 
F nj2 = {A/t)F n -F nA . 



Then, F n ^ is bounded, and the previous proof applies, while F n> i = ^G n , with G n 
smooth, approaching a constant at infinity. Then, 

i[H,F ntl ] = ji[H,A]G n + ji[H,G n ] 

= ±[\p\ + V(x)]G n + j J G n {\)e- iXA ^ J e^/^A/tje-^dsdX. 

The first term is bounded by 2 _n ||G n || on support of E n , and therefore the 
integral over 1 < t < 2 n is bounded by 0(1). The second term is, 

ds 



+ - f I d n (\)e^ t jj e-^[\p\+V(x,a)\da. 

I [ X e s/t d s = - t -e- s / t \ x = \{z~ X/t -^) 
A Jo A A 



and, for f < 1, e~ x l l - 1 ~ -f + 

so the first term contributes J G n (X)e lXA ^ t dX = 0(1), and the second term is 
bounded by cf \XG n (X) \d\/t. 

Furthermore, there is a factor of 2 _n , coming from \p\ + V(x, s) < \p\ + c\p\ 2 . 

Here, V(x, s) = e isA / t V(x) e - isA / t . ■ 
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Section 2. Auxiliary Identities and Inequalities 
Lemma 2.1. Assume H = \p\ + V and \V\ < ^. Then, 



(2.1) 



\p\E In (H)f + VE In (H)f 



< c2~ 



Ei f 

1 n J 



1; 



Proof. For 

(f,E In (\p\ 2 + V 2 )Eij) = (f,E In (H){H 2 -V\p\-\p\V}E In (H)f) 



Let a 



\p\EiJ 



< 2 



-In 



Ei f 

1 n J 



+ 2 



Then, b<(l-6) \p\E In f 



tainty inequality(in 3-dimensions or higher), and so, 

(/, E In (\p\ 2 + V 2 )E In f^) = (a 2 + b 2 - 2ab) + 2ab 



by the uncer- 



It follows that, 



\p\EiJ 



< 2 



-In 



>8< 



Ei f 

1 n J 



\v\EiJ 



+ 2ab. 



Proposition 2.1. // \p\ < H, then (2.1) holds. 

Proof. When \p\ < mH, we have that, by the spectral theorem, 



(2.2) 
Hence, 



1 m 

if - R 



i 



m/ 2 



f 



<m 1 ' 2 



L 2 



1 



\P\ 1/2 ' 



v{x) ¥\ 



< 



V(x)- 2 |p|-l 
r oo \p\ 



< oo, 



V(X)[1 H-U 



= -V(x)^V(x)^- 

tl \p\ 

= —V(x) — -rV(x)r- — 
H r r \p\ 

= -V(x)r-^-r 2 V(x)-^-. 

r H r r \p\ 
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Since, 

we get that, 

V(x) 



111 

r H r 

1 

H 



1 1 



1 1 



r H 1 / 2 H 1 / 2 r 



1 1 



1 1 



1 

bl 



< 



l 



m/ 2 l rV(x] 



1 1 2 ( Al_ 

H l/2 r r V ^ X > r \ p \ 



< m 



1 1 



|p|l/2 r 



i i 



bl 1 / 2 



(r 2 V) 



1 1 



where Cy is a constant, depending on the L°° norm of V, rV, r 2 V. 
We therefore conclude that, 

is also bounded. Finally, we have that, 

|p|£ /n (ff) = \p\^HE In (H) = 0(l)2~ n E In (H). I 

Now we prove some useful identities. 
Lemma 2.2. 

(i) r 2 (-A) = A 2 + L 2 - iA - 3/4, (dimension= 3). 



(ii) e iXA \p\ a e~ iXA = e~ aX 

(Hi) gi^A^ctg—iXA _ ^a\ r a 



\P\ C 



A^Urdr + l). 



H W\ Al r = W\^W\ - f )f = R<V A - t) J = 0(1) 

^ ||^^||<2|^||. 

Proo/. 

(i) r 2 (-A) = r 2 (-a 2 - f a r + = -r 2 a r 2 - 2rd r + L 2 

= -2rd r + 2rd r - rd 2 r + L 2 = -rd 2 r + L 2 , 

A 2 = -\(-2rd r - 3) = (-4rd r rd r - 9 - I2rd r )\ 

= rd r - rdlr - I - 3rd r = -rdlr - 2rd, 
so, 

(ii) d x (e lXA \p\ a e- lXA ^ = 

e iXA \p\ a e~ iXA = e~ Xa \p\ a . 

(iii) As in(ii), but now, i[A,r a 

(iv) Follows from (v). 



_ 9 

r 2 (-A) = -rd 2 r + L 2 = A 2 + § + L 2 - iA 



\=A 2 + L 2 -iA-\. 



e lXA i 



A, \p\ a \e 



,-i\A 



—ae lXA \p\ 



otg—iXA 



ar 



(v) 



\P\ r r 



L 2 



< 2 



r-ip 



L 2 



2 1| by the uncertainty inequality in 3 dimensions. l| 
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Proposition 2.2. 

a) E n {H) = E n {H)E n {\p\)E n {H) + E n (H)8- 1 0(2- n )E n (H). 

b) E n {\p\)E n (H) = E fl (\p\)0(2- n )E n (H), n ^ n. 
For rj > In 2, 

c) E n (\p\)e^ A E n (H) = E n {\p\)0{2- n )E n {H). 

d) E n (\p\)E n (H) = E n (\p\)0(2- n n)E n (H). 

Proof. Part (a). 

E n {H) = -E n {E n {\p\) - E n )E n + E n {H)E n {\p\)E n {H), 
E n (\p\) - E n {H) = jE n (X) (e^IPl - e iXH ) dX 

= J E n (X)e lXH J Q A e- lsH (-i)Ve ts \p\dsdX. 

\V\ < c < x >~ 2 , 

=S> \((j),E n [E n (\p\) - E n (H)} E n ip)\ 

< cj \\< x >" x e~ isH e iXH E n (j)\\ \E n {X)X\\ ||< x > _1 e~ is ^ \E n ij}\\ dsdX 

< C5- 1 \\HE n </>\\ J \XF n (X)\dX \\\p\EM 
<c5- 2 2-^2-^U\\=c5- 2 2- n U\\Ul 
where we used that, 

< x > _1 < c\p\ < c5~ 1 H 

in dimension three or higher, and 
proposition(2.1). 

Notice that, \\pE n if>\\ < 1| || < «y _:l 2- n ||V>||, by Lemma 2.2. 

Part(b). 

\{^E n {\p\)E n {H)^)\ = | (<t>,En(\ P \) [En(\p\) - En(H)] E n (H)lp) \ 

<cS-^HEM\\\p\E n (\pMJ\E n X\d\, 

by the proof of Part (a). The last expression is therefore bounded by, 

C 5- 1 2-||^||2-^||0||2^ = C 5- 1 2-||0||||^||=O(2-"). 
Part(c) follows from Part(b), since, for |A| > ln2: 

E n (\p\)e iXA = E n (\p\)e iXA E n (e- x \p\) 

<E n {\p\)e iXA ^(bl). 

\n\<M 

with n 7^ n, M < oo. 

Part(d) follows from part(b), since the domain of \p\, in the support of E n (\p\), 
is covered by n dyadic intervals from [2 _n , !].■ 
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Proposition 2.3. Assume that H = \p\ + V(x) has no bound states or zero 
energy resonances, and that V(x) vanishes faster than r~ 2 at infinity, and is suffi- 
ciently regular. The dimension is 3. Then, for some < m < oo, 

(2.3) H>m\p\. 

Proof. Since H has no bound states, and V(x) — > at infinity, H > 0. If we now 
make a small perturbation, H — >■ H e = H + eV, then, since H has no zero energy 
resonances, H e has no bound states, for e sufficiently small. Hence H e > 0. 
But then, 

H= Y+~e ( H + (1 + e)v ) + rb) bl = TT7 He + m|p| - mlplu 

Lemma 2.3. For H = \p\ + V , we have that \p\ > 5H , for dimension= 3, 

provided, 

\V\ < c/r. 

Proof. For 6 > 0, 

(2.4) 5H = S\p\ +SV < 5\p\ + 5c\p\ = 6(1 + c)\p\ < \p\, 

where we used that for \V\ < c|a;| _1 , we have |V| < c\p\ in 3 dimensions. 
So, in particular, we have, 

(2.5) E In (H)\p\E In (H) > E In (H)5HE In (H) > 6m£I n E In (H). 
■ 

Proposition 2.4. Suppose, as before, that H = |p|+V(x), the dimension is 3, and 
that V(x) is sufficiently regular, and vanishes faster than r~ 2 at infinity. Suppose, 
moreover, that H has no zero energy resonances, and no zero energy bound states. 
Then, 

(2.6) P C (H)HP C (H) > P c (H)5\p\P c (H), for some 5 > 0. 
Proof. We have that, 

P c (H)(\p\ + v)p c (H)>Q. 

Add a small perturbation eV to H: H e = H + eV = \p\ + (1 + e)V. 
Then, for e sufficiently small, no new bound states are created. 
Hence, 

(2.7) f(H e > -e )(b| + V + eV)f(H e > -e ) > 0, 
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for some eo > 0, and / is a smooth characteristic function of the interval [— eo, oo] 
Then, 

P C (H)HP C (H) = P C (H) (j^He + ei|p|)p c (^), 
would imply that: 

P C (H)HP C (H) > ei P c (#)b|P c (#), 
if we can prove that for some e 2 < ei , 

(2.8) P c (H)H e P c (H) > P c (-e 2 \p\)P c . 
So, want to use (2.7) to prove (2.8). 

P c (H)H e P c (H) = P c (H)f(H £ < -e )H e P c (H) + P c (H)f(H e > -e )H e P c (H) 

(2.9) > P c (H)f(H e < -e )H e P c {H), 

by (2.7). Now, 

P c (H)f(H e < -e ) = P C (H) [f(H e < -e ) - f(H < -e ) 



P c (,e") J f(\)e +iXH i e- lsH eVe isH *ds 



= P c (H)0(eV) J \Xf(X)\d\ 
= P c (H)0(eV)0(e^), 

so, 

P c (H)f(H e < -e )H e P c (H) 

= P c (H)O(eV/e )f(H e < -e )H e (f(H e < -e ) + f(H e < -e j)p c (H) 

= P c (H)O(eV/e )f(He)HJ(H e < -e )P c (H) 

= P c (H)0(eV/eo)f(H e )H e O(eV/e )P c , 

since f(H e < -e )f(H e < -e ) =0. (/ + / = 1, // = /). The last term, can be 
bounded by: 

(2.11) P c (H)e 2 < x >" 2 O eo (H e ) < x >" 2 P C (H) < O eo (l)P c e 2 \p\P c , 

O eo ~ O(l/e )H e , coming from : 



J \Xf(X)\dX = O(l/e ). 



€q is basically the distance of zero to the largest (negative) e.v. of H e , and e is 
arbitrarily small, so t 2 jt\ < 0(e 2 ). Hence, (2.9) - (2.11) imply (2.8). ■ 
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Section 3. Maximal velocity bound 
We begin with estimating, for a > 1, 



(3.1) E n (H)F a \U- > a J E n (H) = E n F a E n 

(3.2) E n F a E n = E n F a (F n (A/t) + F n {A/t))E, 
where, 

(3.3) F n (A/t) = F n (j > R2~ n ) = 1 - F n , 

(3.4) 1< R < a, 
and such that 

(3.5) F n (2- n b)F a (b) = 0. 



Since, by the propagation estimates, 

\\F n (A/t)E n (H)^(t)\\<o(l)\\EM 

with o(l) — > as t — > oo, we need to control E n F a F n E n by a decaying function of 
t, in order to prove the maximal velocity bound, on the energy shell I n . 

Proposition 3.1. 

\\E n F a E n \\ < o(l)\\E n iP\\ + —\\E n iP\\ 

Proof. 

E n F a F n E n = E n (\p\- 1 \p\- 1 A)r' 2 F a r-' 2 F n E n 

= E^-^-^A 2 + cA + Cl ){r- 2 F a )F n E n 
= c 2 E n \p\- l (\p\- 1 Ar- l )(A + c)r- 1 F a F n E n 
(3.6) + csEnlpl-Wpl^r-^r^FaFnEn = B x + B 2 . 
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(3.7) 

where we used that 



B 2 = E n O(2 n )O (l)c 3 - 



E^pl' 1 = 0(2 +n ) and (ppV" 1 = O (l). 



^)F n 



E E 



Since t fF a = 0(l), B 2 = 0(2 n /at), 



B 1 = £ n 0(2")0i(l)c-L (^F a ) F n £, 

+ E n O(2 n )0 1 (l)\ (-F a ) Rt(Rt)~ L AF n E, 



at V r 



+ £ n O(2")0 1 (l)- 



/at 

' (7 F ° 



where 



Therefore, 



Oi(l) = |p| _1 (A - 3/2z)r" 1 = • xr' 1 = £ fp- 



B 1 = E n O{2 n /at)F n E n ^E n O{2 n )0 1 {l) I -F^ -0(2- n )F n E r , 



a 



(3.8) 



0(2 n /at) + ( ^ ) . 



The key to this computation is the repeated use, as we do below, of the following: 
E n = ^M^br 1 ^ 2 + cA + Cl )r" 2 

= c 2 E n \p\- 1 (\p\- 1 Ar- 1 )(A + c')r- 1 +c 3 E n \p\- 1 (\p\- 1 r- 1 )r- 1 

= Enipi^OiWiA + c'y- 1 + caKbr^M-V- 1 )/-- 1 

(3.9) +c 4 E n \p\- 1 (\p\- 1 r- 1 A)r- 1 . 

We apply it again, to the 0(f) term: (F^ (A) = A fc F(A)); 
0{R/a) --EnO^Fi^FWEn 



a 



= -E n \p\- 1 1 (l)(A + c)r- 1 O l (l)Fl^F^E r , 



a, 



(3.10) 



+ -E n \p\- 1 O (l)r- 1 O 1 (l)F^F^E n . 
a 
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Now, the important observation is that, 

r~ x \p\~ l Ar~ l = (r~ 1 \p\~ 1 A)r~ 1 

and 

r~ 1 o 1 {i) = o^iy- 1 - r^o^- 1 )^ 1 = o^iy- 1 + oiiy- 1 

(3.11) 
We derive: 

1 (l)(^ + c')r- 1 1 (l)Fi- 1 ^ 1 ] = 1 (l)(^ + c')(r- 1 p- 1 ^r- 1 )Fi- 1 ] J Fi 1 ], 



Ar-^FW = r-'F^AF^ + cr^F^F^ 

= ±- t Fl-V(Rt)2-FW + ^Fi^FW + l {dxF l-^ F m 
= 0^ t ) + 0{2-^F^FW). 



OF 

Using that Ar~ x p~ x = Oi(l) + O (l), we have, 



-EnOtMFWFWEn = -EM-'Oiiiyr-'lpl-'A + c'y-'Ft^FmEn 



a a 



+ ^E n \p\- 1 1 (l) 2 r- 1 Ft 1 ^AF^E ri 
ct 

2 



(3,2) = (||) 0^?Fl-*m n . 

where we used that, 

E n \p\- X = E n O(2~ n ), r 2 A = A 2 + cA + c\ 

IpI" 1 ^- 1 = 0(1). 
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We now do this computation again, this time, for the term of order (-f ) 2 : 

Doing it k times, we get: 

(3.13) 

£ ) " E n O(l) k Fl-VFWE n + £ E n O (^) (|) ^0(iyP^Ft 1] E n , 
with, 

(3.14) Ft k] = ~Ft k+1] + rd r -Ft k+1] . 

at at 

For k ~ 51nt, and — < 1 sufficiently small, we have that the E n F a F n E Ti 
term is bounded by 0(t _1 2 n ).B 

We can now prove the Theorem on Maximal velocity bound. 
Proof of Maximal Velocity Bound 

Now, 

(iP(t),F 2 a (- t >am)) 

= j2(H- 1/2 m,F2E n H^m)+Q 

n 

= (#- 1/2 m F 2 (F n (A/t) + F n (A/t))E n H 1 ^ 2 tp(t)^ + Q 

n 

(3.15) 

= J2( H ~ 1/2 ^ FlF n E n H^) + ff" ^ V, F 2 F n E n H l / 2 ip) + Q. 



n 

< WFaH-V^W ^ < n >- 1 / 2 - £ / 2 < n >V2+ £ /2 

n 

< cWFaH-^W ^ < n >1+£ WFuEnH 1 ' 2 ^ 

< cWFaH-^W (j2o(l)\\H 1/2 E n < X >V2 ^||2 < n >1+A 



(3.16) 

< C ||F a if- 1 /2 V ,|| (l)|| <x >V2 
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J2(H- 1/2 ^FtF n E n H^) 



< 



+ 



^2(F a H~ 1/2 ^, F n F a E n H 1/2 ip) 

n 



(3.17) 



<c|ktf-VV||^ 112^2^ 



1R 

t a 



It follows that, 
(3.18) 

(i>(t),F^(t)) <c FaH- 1 ' 2 ^ \\nH\^H- 1 ' 2 <x> 1 ' 2 Vo o(l) + Q, 



(3.19) 



Q = {H-^m,[H 1/2 ,Fi]m) 



To control Q, (3.19), we need to commute fractional powers of H. To this end we 
use that: 

"oo \a — l 

H a = c n I - -HdX, 



o A + ff 



and estimate, 



A 



a-l 



H 



H + X 



,F a 



= X a ~ x H 
-H 



H + X 



,F a 



+ \H,F a ] 



X 



a-l 



H + X 



t a H + X 



H + X \t a J H + X 
(3.20) dX® = 0(1)0 0(H-«). 

Therefore, using (3.20) with a = 1/2, we have that 

\q\ = \(H-v*mMW(])H- i/2 m)\ < -jH-^mw 2 <- t \\< x > 1/2 



End of Proof. 
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